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FLOWN ET: A COMPUTES PROGRAM FOR CALCULATING SECONDARY 

FLOW CONDITIONS IN A NETWORK OF TUFBOMACH INERT 


James R. Rose 
Lewis Research Center 


ABSTRACT 


A coaputer prograa that calculates the secondary flow 
conditions in a network of turboaachine components is 
described. The types of flow components that can be treated 
are face seals, narrow slots, and pipes. The program is 
written in both structured FORTRAN (SFTRAN) and FORTRAN IV. 
The program must be run in an interactive (conversational) 
mode. 
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ttOWHET: 


A COMPUTER PROGRAM FOP CALCULATING SECONDARY 


FLOW CONDITIONS IN A NETWORK OF TUPBOBACHINFPY 


James P. Rost 
Lewis Psstarch Center 

SUMMARY 




The computer program dascribed in this report calculates the 
secondary flow conditions in a network of turbomachine 
components. The types of components that can be treated are 
face seals, narrow slots, and pipes. The program can treat 
networks containing up to fifty flow components and 
twenty-five unknown network pressures. 

The program requires the network parameters, the flow 
component parameters, the reservoir conditions, and the gas 
properties as input. It will then calculate all unknown 
pressures and the mass flow rate in each flow component in 
the network. 

The network solution is obtained using a least squares 
technique. The mathematical models of the least squares 
technique and the fluid flow are not contained in this 
report. They are dascribed in publications which are 
referenced in this report. 

A portion of the computer program is written in structured 
FORTRAN (SFTRAN) . The rest of the program is written in 
FORTRAN IV. 

The program must be run in an interactive (conversational) 
mode. A sample problem is included. 


INTRODUCTION 


The computer program described in this report calculates the 
secondary flow conditions in a network of turbomachin* 
components. The types of components that may be considered 
are shaft face seals, narrow slots, and pipes. The 
mathematical model for the flow conditions is described in 
reference 1. The user may add other flow models to the 
program. The way in which this is done is described in this 
report. 
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REPRODUCIBILITY OP THE 
ORIGINAL PAGE IF 1* H »H 

Tht approach used to obtain the network solution Ip an 
iterative procedure which relies on a least squares search 
technique. The aatheaatical models employed by the least 
squares search technique are presented in this report and in 
reference 2 . 

FLOWNET is comprised of an executive program, a flow model 
subroutine, a least squares search subroutine, and an 
interpolation routine. The flow model subroutine is a 
modified version of the computer program described In 
reference 1. The executive program is written in structured 
FORTRAN (SFTRAN). The rest of the program is written in 
FORTRAN IV. 

FLOWNET has been run on an IBM 360/67 time sharing system at 
the Lewis Research Center (ref. 5). This report is intended 
to serve as a user’s guide for the computer program FLOWNET. 


SAMPLE FLOW NETWORK 

Figure 1 shows a relatively simple flow network. There are 
two flow paths in the network. The two paths are connected 
in series between two pressure reservoirs (nodes). The 
assumed mass flow direction is indicated by the arrows and, 
of course, this will be the direction of the mass flow if 
P 2 > Pi > P 3 . An integer has been assigned to each node. 
This designates the inlet *iid outlet node numbers of both 
flow paths. 

As an example of the application of FLOWNET to this network 
assume that the pressures at nodes 2 and 3 are known 
(fixed). Assume, also, that no mass is lost from the 
network. This requires that the mass flow into each node be 
equal to the mass flow out of that node. In this example, 
the problem is to solve for the pressure Pi at node 1. The 
program FLOWNET solves for the particular pressure at node 1 
which causes the net mass flow at node 1 to be zero. In 
obtaining this solution the total mass flow in the system is 
also learned. 


ASSUMPTIONS 

To find the solution for a given flow network, several 
assumptions are made: 

1. Each node (a point in a flow network where one or more 
flow paths have an inlet or outlet connection) can be 
treated as a large reservoir in which flow velocities 
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are negligible. There is, thus, no ambiguity In th» 
meanings of node temperature and pressure. 

2. K unique, steady-state solution exists. 

3. The fluid temperature at each node is known and there is 
no heat transfer in the network. 

These assumptions allow the problem to be reduced to the 
solution of a set of simultaneous equations, which express 
the conservation of mass flow at each node in the network. 
The method used in FLOWNET to solve this set of equations is 
described in the next section. 


» 


t. 


METHOD OF SOLUTION 


The method used in FLOWNET to find network solutions is 
described in this section. For each flow path in the 

network let the mass flow be given by the function 
»k ' T )• The subscript i refers to the flow path’s 
inlet node and j refers to its outlet node. If all the 
pressures P, and temperatures T, are specified then the 
functions tii.i may be obtained by use of a flow module (a 
subroutine which calculates the mass flow for a specific 
type of flow path) . ht any node K, which is neither a 
source nor a sink, the mass flow equation is obtained by 
adding all mass flow in flow paths having the Kth node as 
the outlet node and subtracting all mass flows in flow paths 
having the Kth node as the inlet node, that is. 


R 


■ 

K 


E 

all i 



E 



( 1 ) 


where r k is called the residual at the Kth node. 
Conservation of mass flow requires that R* = 0 for all K. 
FLOWNET obtains the network solution by solving for the set 
of pressures which results in =0 for all K. 

The strategy used in FLOWNET to find this network solution 
is the least-squares method. In this method the sum of the 
squares residuals. 


all K 

is calculated. Then a systematic search is made in 
parameter space, here the unknown node pressures, for a 
minimum in the sum of the squares residuals, that is. 



R = 


E 

all K 



a minimum 


( 2 ) 
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The first set of mass floe functions generated and used to 
satisfy equation (2) yields only the first approximation to 
the correct values of the unknown pressures in the flow 
network* This first set of mass flow functions is generated 
using user supplied estimates of the unknown pressures. 
Thus, they are not, in general, the true functions for the 
balanced flow network. Thus, iterations must be performed 
until the values of the unknown pressures found by the 
least-squares search routine are nearly equal to ♦■hos» 
values used to generate the mass flow functions. The 
iterations continue until this self consistent solution is 
found. 


node! Function for Hass Flow 


The least squares search routine used by FLOWNET requires 
the partial derivatives of the mass flow functions. The 
mass flow functions are obtained from rather complicated 
algorithms and, thus, cannot be differentiated analytically. 
The numerical differentiation of the mass flow functions is 
simplified by introducing a simple model function 
F i j ( P i/P j) defined as 


4 






( 3 ) 


where the pressure ratio P r is equal to P^ /P* . with this 
model function the partial derivatives of tne mass flow 
functions can be written in terms of ordinary derivatives of 
the model function with respect to the pressure ratio P r . 


GENERAL COMMENTS ON THE PROGRAM 

The executive program is written in Structured FORTRAN 
(SFTRAN) reference 3. An SFTFAN pre-compiler is used to 
generate a FORTRAN IV version of the executive program. The 
rest of the program is written in FORTRAN IV. The entire 
program, FLOWNET, has been run on an IBM 360/67 time sharinq 
system at the Lewis Research Center. Flow networks 
containing several types of flow paths can be treated by 
FLOWNET. Subroutine QUASC (ref. 1) is used to treat flow 
across face seals, through narrow slots, and through pipes. 
The user may also supply his own flow module, provided that 
the module be made to conform to certain interface 
requirements that are delineated later in this report. 


5 



Flow of Program logic 


The general flow of th* program logic is shown in fiqur* 2. 
The executive program reads the data required to define th* 
flow network. This data is then used by the executive 
prograe to model the network for the rest of the program. 
The executive program then makes a series of calls to the 
appropriate flow module to calculate the mass flow functions 
Wi rj (P i»P j#T) for each flow path in the flow network. To 
calculate^ these mass flow functions, initial estimates of 
the unknown pressures must be used. The user must supply 
these estimates. When the calculation of all mass flow 
functions is finished the executive program calls the 
non-linear, least-squares search routine. This routine uses 
the mass flow functions and their derivatives to solve for 
that set of unknown pressures which minimizes the sum of the 
squares residuals (see equation (2) ) . This set of pressures, 
however, is only the first correction to the user’s 
estimates of the unknown pressures. To obtain a self 
consistent solution to the network an iterative scheme is 
required. This iterative scheme is described in the next 
section. 


Iterative Schemes Used By FLOWNBT 

In the previous section, the first pass through an iterative 
loop was described. The iterations continue as follows, see 
figure 2. The first set of unknown pressures calculated by 
the least-squares search routine are used to calculate a new 
set of mass flow functions. Then the least-squares search 
routine is called again to calculate a set of pressures that 
satisfies equation (2) for this second set of mass flow 
functions. This process continues until a particular 
built-in convergence criterion is satisfied, namely. 


<i) _ p (i+l) 
K K 


< 


<AP) 


max 


(4) 


where, 

the superscript refers to the iteration number, the 
subscript refers to the node number and (&P), max is program 
input. When relationship (4) is satisfied for each unknown 
network pressure the solution is considered to be converoed. 
Note in figure 2 that the least-squares search routine is 
iterative also. This is because of the non-linear nature of 
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th« unknowns. 



MISCELLANEOUS FEATURES OF FLOWNET 

Several features have been built into FLOW if** which 

(1) eauaa tha nass flow functions to ba iota accurate in 

tha vicinity of tha solution aftar aach outar iteration 

(2) provide automatic treatment of flow reversal, and 

(3) parsit tha usar to traat ‘.char typas of flow components 

by addin? his own flow aodula to tha program. 


Solution Accuracy 

In general, the converged solution will ba accurate if all 
final nass flow functions are accurate in tha vicinity of 
tha solution. Tha usar cannot insure this accuracy with 
input data (or any other feasible way) because this would 
require having a good knowledge of tha converged solution, a 
priori. Tha first set of nass flow functions are generated 
using estimates of the unknown network pressures and a set 
of input pressure ratios. These pressure ratios nay yield 
nass flow functions which are not accurate in the vicinity 
of the converged solution. This difficulty is overcome as 
follows. After each inner iteration, updated values of the 
unknown network pressures are available. These pressures 
are used to calculate a new pressure ratio for each flow 
path in tha network. Then before generating the mass flow 
functions for the next outer iteration this new pressure 
ratio is placed into the existing pressure ratio array and 
an existing pressure ratio value is discarded. hus, after 
each outer iteration the mass flow functions are more 
accurate in the vicinity of the solution. 


Flow Reversal 


The program is also designed to handle so-called flow 
reversal. This situation occurs when the user inadvertently 
reverses the definition of inlet and outlet nodes in the 
input data so that P A _ < ^ ut . This can also occur if the 
user expects a particular node, where the pressure is 
unknown, to he an inlet node and the calculation (converged 
solution) shows it to be an outlet node. Again P in < P 0 ut 
is obtained. The program handles these situations by 
interchanging the definitions of the inlet and outlet nodes 
cl the path or paths in question. The user is informed that 
these definitions have been changed. 


User Supplied Flow Module 

FLOWNET has a provision for the user to supply an additional 
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flow nodule. Thu user nay want to consider flow conponen*« 
having flow characteristics that art not described by ouasc 
( caf. 1). A subrontlna that describes the desired flow 
charaetariat lea aay than ba added by tha uaar. Thin user 
supplied flow nodule auat ba aada to lntarfaca with FIOWNPT. 
Thla lntarfaclng la achiavad at followa. 

(1) Tha flow aodala auat ba eallad subpoutine ADDPlo. 

(2) Tha f 1 o* f aodula auat not hava an arguaant vactor. 

(3) Tha following labeled COHMON statement nust appear In 
tha aubroutina. 

COHHON/NODULE/ONLYI N, NDONLY , PINLET, PEATIO,TO, INPUT, PLOW, TP A TH 

Tha valuaa of all quantities In cohhon/hoddlb/, except plow, 
ara aatabliahad by tha axaeutiva prograa, Tha valuaa of 
thaaa guantitiaa dataralna: (1) Tha purpoaa of tha call to 

tha flow aodula by tha axaeutiva prograa (thare are three 
purposes) and (2) the values of tha inlet pressure PINLET, 
tha pressure ratio PRA ¥ 'TJ, and the teaperature T0 r to be 
used by tha flow aodula to calculate the aass flow plow. 

On a given call to tha flow aodula, tha purpose of tha call 
is deterained by the values of tha integer INPUT, and the 
logical ONLYIN. Whan INPUT equals one and ONLYIN is .TPUE. 
tha flow aodula aust read the flow path input data and write 
this data into tha high spaed printer output dataset. No 
calculations need ba dona on this call. When INPUT equals 
one and ONLYIN is .FALSE, the flow nodule must read the flow 
path input data and calculate the aass flow, when input 
equals three and ONLYIN is .FALSE, the flow nodule 

calculates the nass flow. The flow path input data nust not 
ba read on this call. The logical NDONLY always has the 
value .TRUE.. The use^s flow nodule nay do more than 
calculate aass flow. It nay, for exanple, calculate radial 
profiles. NDONLY nay be used to skip unnecessary 
calculations and/or output. An exanple of how the 
quantities in COHHON/NODULE/ are used in a flow nodule i« 
shown below. 
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SAHPLE FLOW HODOLE 


SUBROUTINE ADDFLO 
LOGICAL ONLYIN, HDONLY 

COHNON/HOD OLE/O WLYI B , HDONLY , PIN LET, PP AT 10, TO# INPUT, FLOW.IPATfl 

t 


IF (INP0T.BQ.1) GO TO 5 
IF (INP0T.BQ.3) GO TO 10 
BEAD THE FLO* PATH INPOT DATA 
IF (.NOT. ONLYIN) GO TO 20 
WRITE THE FLOW PATH INPOT DATA 
GO TO 25 
CONTINOE 


t 

PO * PINLET 
P3 « PO/PRATIO 
TEHP - TO 


HASS FLOW - F (P0,P3,TEHP) 
FLOW * HASS FLOW 
IF (HDONLY) GO TO 25 


(Other calculations and/or output not desired for flow 
network usage) 


25 CONTINOE 
RETURN 
END 


5 

10 

15 

20 


10 


IN POT OITA 


The Input data required by PLONNET ara of two genaral typaa. 
Tha firat typa la the flow network data. Thin data la read 
by tha executive prograa. other data tea ralatad to 
apaciflc flow patha In tha natwork. Thia flow path data la 
raid by subroutine QUASC. 


Natwork Data 

Tha natwork input data la raad by tha aain prograa. This 
data auat ba raad froa unit 4. Tha first card In tha 
natwork data la a titla card. Tha titla card la uaad to 
idantlfy tha particular natwork to ba traatad. card colunna 
1 to 72 of ona card nay ba uaad and thla card la raad by 
format (18*0). Tha rant of tha natwork input data la in 
naaaliat foriat. Thara ara thraa naaallata with tha naaaa 
PATHSP, NODESP , and PARAN. Tha carda following tha titla 
card have tha data in N A H EL 1 5 T/PA TH SPA Thia data ia listad 
in TABLE I. Tha naxt carda have tha data in 
NAHELIST/NODESP/. Thia data ia liatad in TABLE II. Tha 
last cards in tha natwork input dataaat hava tha data in 
NAN ELIST/P ARAN/. This data ia liatad in TABLE III. 


Plow Path Data 

C'ha flow path data auat ba raad froa unit 5. This data is 
raad by subroutina QUASC. Tha first card required by 
subroutina QUASC ia a titla card. Tha titla idantifiaa tha 
data for tha given flow path. Tha titla card usas columns 1 
to 72 of ona card and ia raad by foriat (18A4). Tha naxt 
cards required by subroutina QUASC contain tha paraiatars in 
NANELIST/SDATA/. Thasa parameters ara listad in TABLE iv. 
Sons of tha paraiatars in NANELIitT/SDATA/ diffar from thosa 
in tha varsion of QUASC dascribad in rafaranca 1. 

Also, some of thasa quantitias must always ba the same value 
in tha modified varsion of QUASC uaad in this report, whan 
this is tha case., TABLE IV shows tha value that must he 
used. Tha naxt card contains tha integer NJ, which is the 
number of film thicknesses to be considered for a given flow 
path. When solving flow network problems, only one film 
thickness per flow path can be considered. Thus the value 
of NJ is always 1 and is raad by format (T3) . The last card 
in the flow path data for subroutina QUASC has the value of 
the film thickness and is raad by format (F12.6) . THp 
NAHELIST/PDATA/ used in program QUASC (ref. 1) is not used 
in subroutina QUASC. 
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A sst of the lata describe! in this section mast be included 
in the floe path dataset for each flow path in the flow 
network. 


SAMPLE NETWORK DATA 

Figure 3 shows a relatively simple flow network. The arrows 
indicate the direction of the mass flow in the network. The 
circled numbers are the Integers arbitrarily assigned *o 
each flow path. Mote that there is only one node where the 
pressure is unknown and this must be numbered node 1. i>j is 
the pressure at the Jth node. is the total temperature 
of the gas. The first estimate of the pressure at node 1 
may be obtained by assuming a linear pressure distribution 
across paths 1 and 3, if a better estimate is not available. 

Table V is the set of network data which would be used to 
describe the network of fig. 3. 

PUNNING THE PBOGFAM 

In this -action it is assumed that the reader is familiar 
with the IBM 360 time sharing system (ref. 5). To initiate 
running FLOWNET it is necessary to datadef the required 
input and output datasets and the joblibrary that contains 
the program object modules. The flow network input must be 
read from logical unit tt. The flow module input must be 

read from logical unit 5. Logical units 8 and 9 are used 
for keyboard input and output, respectively, to interact 
with the program. Logical unit 6 is used for program output 
obtained from the high speed printer, A typical set of data 
definitions are: 

ddef ft04f001, vs, network. input 

ddef ft05f001 ,vs, flow. path. input 

ddef ft06f 001, vs, network, solution 

ddef library, vp, object. f 1 owne t, option =joblib 

Logical units 8 and 9 are defaulted to the keyboard, After 
the above data definitions have been performed, program 
execution is initiated by simply entering the following at 
the keyboard, 

USER: CALL MAIN (name of the main program object module) 

This command causes the first set of mass flow functions to 
be calculated, i.e., the first outer iteration is initiated. 
One mass flow function is calculated for each flow path in 
the flow network. These functions are calculated using the 
known pressures and temperatures in the flow network, the 
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unr' j«tinat«s of tha unknown pressures, and the input 
values of the pressure ratios. After these functions are 
calculated the following prompt appears at the keyboard, 

HAIIf: OUTPUT PLOW FUNCTIONS POP THIS ITERATION? 

The user then responds, 

USER: y (or n) 

to obtain (or suppress) the output of the nass flow 
functions for this iteration on the high speed printer. 
After this, the non-linear least squares search routine, 
NOWLIN is called by the main prograa. This initiates the 
inner iterations. The user is proapted to interact with 
NONLIN as follows, 

NOWLIN: ENTER (1,2) POR (CONVERSATION, BATCH) 

USER*. 1 

Note that NONLIN is a general purpose prograa. However, 
when used with FLOWNET it can be used in a conversational 
(interactive) mode only. Thus, the user always replies with 
the number one to this prompt. 

NONLIN: NEW STARTING VALUES ONLY? 

USER: no (or yes) 

The user must respond "no” the first time this prompt is 
given. A "yes" response indicates a desire for restart or 
failure-recovery procedures. In this case NONLIN would 
immediately prompt for a new set of non-linear parameters; 
here the unknown pressures in the flow network. The "no" 
response, however, causes the following prompt. 

NONLIN: ENTER NF, NR, NQ, NN, OR QUIT 

(For a problem having only non-linear parameters, the 
unknown network pressures in this case, NF and NN are the 
only quantities that have non-zero values. NN is the number 
of data points to be fit. NF is the number of non-linear 
parameters being used to fit the data. Thus NN and NF have 
the same values. This value is equal to the number of 
unknown pressures in the flow network.) 

USER: Sdata nn=2, nf=2 Send 

{This response is for a flow network with two unknown 
pressures. Note the use of NAMELIST/DATA/. Had the 
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response been (Bdata quit*, true. Send) , the inner iteration 
loop would have boon toroinatod and a now oot of oaos flow 
functiono would havo boon calculated.) 

Tho prompts and typical responoes proceed as follows. 

NONLIN: ENTER INITIAL NON-LINEAR PARAHETBRS 

USER: Bdata phi(1) «200.0, phi (2) *25.0 Bond 

(Tho phi's are tho current ostiiates for the unknown 
pressures. Again, note tho use of NAMELIST/DATA/) 

NONLIN: ENTER (1,2, 3,4) POR (GRAD, LIN, SPIRAL, CONP2) 

USER: 3 (or 1 or 2 or 4) 

(GRAD, LIN, SPIRAL, and COMP2 are four different search 

techniques. The search tohniques are described in Appendix 
A and reference 2. The user say use any of then, however, 
the SPIRAL method has generally been the most successful on 
the problems solved with FLONNET at LeRC.) 

At this nonlin begins search procedures according to 

the method selected. Each tine NONLIN determines a new set 
of PHI values (unknown pressures) an interchange takes place 
of which the following is typical: 

NONLIN: R * 1.40524 B-3 

PHI = 2.22054 E 02 1.64304E 01 

ANSWERS GOOD ENOUGH? 

USER: no (or yes) 

NONLIN: ENTER (1,2, 3, 4) FOR (GRAD, LIN, SPTR AL , C0MP2) 

USSR: 2 (or 1 or 3 or 4) 

This interchange continues until the user's response is 
"yes", in which case the final values of R and PHI are 
displayed as follows, 

FINAL RESULTS 

R = 1.876544 E-09 

PHI = 2.19386 E 01 
1.62569 E 01 
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where R is the value of the sum of the squares residual 
(eqn. 2) and the PHI'S are the related valueu of the unknown 
network pressures. (These are not necessarily the final 
results for the converged network solution. They are the 
results at the end of the current inner iteration.) The 
user is then pronpted with, 

NONLIN: NEW STARTING VALDES ONLY? 

user: no 

NONLINs ENTER NF, NP, HQ, NN, OR QUIT 
USER: Gdata quit * .true. Send 

At this point, the inner iteration is stopped and proqran 
control returns to the executive progran. A check is sade 
there to see if the current values of the unknown pressures 
are such that they may oe considered the converged solution. 
If so, this solution is printed out and the program cones to 
a terminated stop. Alternatively, if the current values of 
the unknown pressures do not satisfy the convergence 
criterion, equation (1) , then the next outer iteration is 
begun by calculating the next set of mass flow functions. 
These functions are generated using the set of values of 
unknown pressures just calculated by NONLIN (the last inner 
iteration). The procedure just described continues until 
the convergence criterion is satisfied. 


OUTPUT 

Because FLOWNET is run in an interactive mode there is 
program output at the computer terminal in addition to the 
high speed printer output. The computer terminal output and 
the high speed printer output are described in the next two 
sections. 


Computer Terminal Output 

The purpose of the computer terminal output is to give the 
user information about the status of the solution at the end 
of each outer iteration. The first output at the terminal 
tells the user how closely the mass flow conservation 
equations are satisfied for the current set of mass flow 
functions at each node where the pressure is unknown. This 
output consists of the following information at each node 
where the pressure is unknown. 

1) HD (NET) - The net mass flow (=0. when the mass 
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conservation agnation is satisfied exactly.) 

2) HD (HIM) - The magnitude of the smallest aass flow 
component. 

3) HD (MIT) /MD (HIM) - The scaled net mass flow. 

The values of the scaled net mass flows are particularly 
helpful. They are the best measure of how closely the mass 
flow conservation equations are satisfied. The sum of the 
squares residual R, and the values of the net mass flow can 
actually be misleading if used to decide how good th* 
solution is because they are dimensional quantities. 

The next output at the terminal is the scaled values of th* 
pressure changes from the beginning of one outer iteration 
to the end of the corresponding inner iteration. This 
quantitiy is called DBLTAP/p and is equal to the left hand 
side of equation (6). Its value at each node where the 
pressure is unknown is displayed at the terminal. 


High Speed Printer Output 

The high speed printer output, hereafter called the computer 
output, consists of tha input data, calculated parameters, 
intermediate network solutions, and the converoed network 
solution. 

The first page of computer output is the input referred to 
previously as the network input (See TABLBS I, II, and tit.) 
The next output is tha input and calculated parameters for 
the flow paths. There is one page of flow path output for 
each flow path in the flow network. This output is 
described in detail in reference 1. 

The computer output immediately following the flow path 
output is printed each iteration. This is the output for 
the intermediate network solutions. It consists of, 

(1) Mass flow function information and; 

(2) Mass flow balance information. 

The mass flow function information is printed each iteration 
(if the user asks for it) . The first line of this output is 
the path number, the inlet pressure used to calculate the 
path's current mass flow functions P(INLFT), and the path's 
current outlet pressure P(OUTLET). Following this line of 
output are the current values of the mass flow function 
MASS FLOW, and tha values of the pressure ratio used to 
generate the mass flow functions PPATIO. This mass tlow 
function information is output for each flow path in the 
flow network. 
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The Baas flow balance infotaation appears next in the 
cooputac output* This output la printed at tha coapletion 
of aach innai Iteration and is, also, part of tha coaputer 
tarainal output. This output consists of tha nat aass flow 
at aach nods HD (NET), tha saallast aass flow coaponant at 
tha noda HD (HIM) * and tha acalad aass flow at tha noda 
HD (NET) /HD (Hilt) . 

Tha final coaputar output is tha converged solution for tha 
flow network. This consists of tha following for aach flow 
path in tha flow network; tha path nuabar, tha inlat noda 
nuabar INNODE, tha outlat noda nuabar OOTNODE, tha inlat 
prassura P (INLET) , tha outlat prassura P (OUTLET) , and tha 
aass flow through the flow path HASS PLOW. A saaple problea 
is given in Appendix B. 
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APPENDIX A 


PROGRAM DESCRIPTION 

This appendix contains a description of each of the four 
progress that conprise FLOW NET. A brief, general 
description is given for the flow sodule ODASC (ref. 1) , the 
non-linear, least sguares search routine NONLIN, and the 
interpolation routine INTERPOL used by FLOWNET. A sore 
detailed description of the executive proorras is given. 


DESCRIPTION OF THE EXEC DTI VE PROGRAM 

The executive program is written in Structured FORTRAN 
(SFTRAN) . It consists of a sain program and 7 subprograms. 
This section contains a functional description of each 
subpart of the executive program. 

Main Program 

The main program first reads in the flow network data. It 
then uses this data to construct arrays that provide a 
program description of the network geometry. After this, 
the main program provides the controlling logic for the 
inner and outer iteration loops and checks for convergence 
of the network solution at the end of each outer iteration. 

Subroutine MDFONC 

This subroutine is called by the main program. It contains 
the controlling logic for calculating the mass flow 
functions for all flow paths in the flow network. After the 
mass flow functions have been calculated the user is 
prompted to enter Y(YES) or N (NO) to obtain or suppress the 
print out of the mass flow functions as part of the high 
speed printer output. The subroutine is called once per 
outer iteration. 

Subroutine FPRDPR 

This subroutine is called by the main program. It 
calculates an array of ordinary derivatives of the auxiliary 
functions F if j (P r ) for each flow path in the flow network. 
The auxiliary function Fi,j(Pr) is defined by the model mass 
flow function, 


M 


P . 
1 


i* j 




(A 1) 
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These derivatives art obtained by a numerical 
differentiation scheme contained in the interpolation 
routine IHTEEPOL. The ordinary derivatives of fi,j< p r) are 
calculated at each value of tho pressure ratio currently in 
the pressure ratio array. The partial derivatives of the 
■ass flow functions are required by the non-linear, least 
squares search routine NONLIH. They are obtained from the 
■odel aass flow function and the ordinary derivatives of the 
auxiliary function. This subroutine is called once per 
outer iteration. 


Subroutine YVALUB 

This subroutine is called by the non-linear, least squares 
search routine. It Is used as a leans of connunication 
between the executive progran and the general purpose 

least-squares search routine NONLIN. The subroutine passes 
the following infonation to the search routine for all 
non~fixed nodes in the network: 

(1) The input values of the net mass flow. 

(2) The calculated value of the net mass flow for values of 

the unknown pressures selected by the search routine. 
These values are obtained by interpolation on the 

arrays of mass flow functions. 

(3) The partial derivatives of the mass flow functions for 
values of the unknown pressures selected by the search 
routine. These values are obtained by interpolation on 
the arrays of the partial derivatives of the mass flow 
functions. 


Subroutine close* 

This subroutine is called by the main program. It is called 
after each outer iteration. At the end of each outer 
iteration new iterates for the unknown pressures are 
available. These new iterates are used in this subroutine 
to calculate updated pressure ratios for all flow paths in 
the network. These updated pressure ratios are then 
inserted into the existing pressure ratio arrays. An 
existing pressure ratio value is discarded from the array. 
The pressure ratio values for each flow path remain 
monotonic increasing in the array. 

subroutine flows 

This subroutine is called by SUBROUTINE YVALUE. It computes 
a component of mass flow at a particular node each time it 
is called, rt also assigns the correct algebraic sign to 
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the mss flow. Tha magnitude of the flow in a given flow 
path is determined by intarpolation in the corresponding 
Bass flow function array. Tha aubroutina is called K times 
for aach non~fixed node, whara K is tha nuibar of flow paths 
with inlat or outlat at tha node. Thasa nass flow 
cosponants ara used by SUBROUTINE YYALUE to calculata tha 
nat sass flow at aach node for tha sat of pressures 
currently being used by tha least-squares search routine 
NONLIN. 


Subroutine PDBRIY 

This subroutine is called by SUBROUTINE YYALUE. At a given 
non-fixed node it calculates all required partial 
derivatives of tha mass flow functions. One partial 
derivative is calculated per call to SUBROUTINE PDBEIV. The 
algebraic sun of thasa partial derivatives is the desired 
partial derivative of the net nass flow function at a given 
node. 


subroutine BALANS 

This subroutine is called by the main program. The 
subroutine calculates and prints mass flow balance 
information at each node where the pressure is unknown. For 
each such node the scaled and unsealed net mass flow is 
calculated and printed out. The scaled net nass flow at a 
node is egual to the net mass flow divided by the smallest 
coaponent of the mass flow at the node. 


DESCRIPTION OF SUBROUTINE OUASC 

A brief, general description of the flow module QUASC is 
given in this section. Reference 1 contains a complete 
description of the program. 

Subroutine QUASC is a modified version of the computer 
program described in reference 1. It is called by the 
executive program. It calculates the mass flow of a viscous 
fluid through narrow passages with parallel walls. However, 
if one assumes that the equivalent hydraulic diameter of a 
narrow slot equals twice the distance between walls, then 
QUASC can also be used to calculate flow in other types of 
ducts with constant cross-sectional area. For example, a 
pipe (circular cross-section) with radius r may be assumed 
equivalent to a slot with a gap equal to r and a width egual 
to rrr (so as to have the same flow nrea) . 

The mathematical model used to calculate the mass flow in 
subroutine QUASC is developed in reference 4 and summarized 
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ORIGINAL PAGE IS POOR 



in reference 1. Tht analysis includes fluid inertia, 
viscous friction, and sntrnncs lossss. Ths model is valid 
for both completely subsonic flow and choksd flow, it is 
also valid for both lasinar and turbulent flow regimes. 

It is assuned that tha flow in tha seal leakage flow region 
behaves as a constant-area adiabatic flow with friction, k 
guasi-one-diaensional approximation is nade wherein it is 
assused that the flow properties can be described in terns 
of their cross-sectional averages. 

The following assumptions were nade in the analysis: 

1. The area expansion due to radius Increase is neglected. 

2. The flow is adiabatic. 

3. No shaft work is done on or by the system. 

4. No potential energy difference is present such as caused 
by elevation differences. 

5. The fluid behaves as a perfect gas. 

€. The sealing surfaces are parallel. 

With these assumptions, the flow is commonly known as Panno 
line flow. 

The set of equations in the mathematical model car.not be 
solved explicitly. Thus, iterative procedures are used. 

There are three types of flow considered: (1) critical 

flow, when the exit Nach number is 1 and the exit pressure 
equals the ambient pressure; (2) supercritical flow, when 
the exit Nach number is 1 and the exit pressure is greater 
than the ambient pressure; and (3) subcritical flow, when 
the exit Nach number is less than 1 and the exit pressure is 
equal to the ambient pressure (flow is entirely subsonic) . 

Because the boundary conditions on the mathematical model 
are slightly different for each type of flow, the solution 
of the equations is slightly different for each type. 


DESCRIPTION OF NONLIN 

NONLIN is a general purpose, least squares search routine. 
In the program FLOWNET, this routine is called by the 
executive program to perform the inner iterations (see 
figure 2) . The user is first prompted for input data. 
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Including initial valuta for tha if? non-linear oarameters. 
That* non-linear piraaatars art callad PHI in N3NLTN. ?h» 
uaar it than invited to chooaa one of four poaalbla search 
techniques. These aethods ate called GRAD for the qradlent 
method, LIN for the linear method, spiral, (reference 2) a 
coaproaise aethod and C0HP2 for the second coaproaise 
net hod. 

NONLIN firat deteraines the lowest possible sue of squares 
residual R using the initial paraaeters PHI, and then 
searches for a still lower value by aodifying the PHI values 
according to the search technique selected. 

If the sethod GRAD is selected, NONLIN calculates a vector 
which points in the direction of steepest descent froa the 
current point P in non-linear paraaeter space. An initial 
step is taken in this direction: if P decreases, another 
step is taken: if R increases, the distance is halved and R 
calculated again. This process continues until, in the 
first case, R begins to increase, or, in the second case, a 
value is found which is lower than at the initial point P. 

If the aethod LIN is selected, NONLIN calculates a vector 
which, when added to P, gives the point L which would yield 
a nininun R if the linear approximation were valid. Tf the 
value is higher than at the point P, however, backward steps 
along the line LP are taken until either a lower value is 
found or the method is judged to be fruitless, in which case 
the method GRAD is invoked. 

Two compromise methods are also available. If SPIRAL is 
selected, points are axamined along a curve in parameter 
space which begins at L and swings out so that it approaches 
P along the direction of steepest descent. (At the same 
time NONLIN examines a second curve, which is the mirror 
image of the first with respect to a line bisecting the 
angle between the steepest descent vector and the linear 
vector.. The second curve is examined until the two 
intersect.) This aethod is useful when the path of the 
minima leading from the point p to the actual minimum point 
curves and crosses the line PL. If the first attempt fails, 
the starting point (measured from P) is taken to be half the 
original distance, and so on until four spirals are 
examined. If all four attempts fail, the method GPAD is 
invoked. 

The second compromise method, C0MP2, begins by finding a 
lower value of R along the line of steepest descent. Then a 
second lower value is sought, this time along the line PL. 
If both attempts are successful, the line joining the two 
points is examined to determine a still lower value; if not. 
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isprovesent along the line of steepest descent i* sought by 
■tans of REFINE. Jill four methods eventually sake us* of 
REFINE, which Is designed to examine th* sus of squares 
along a glv*n llna PP* to detereln* a "best" value, In th* 
casa of GRAD, LIN, and SPIRAL, th* point p Is th* current 
point and th* point P* has b**n d*t*rsin*d by th* s*l*ct*d 
method; in th* cas* of C0HP2, P has th* hlgh*r and P* th* 
low*r valu* of i for th* two points found. First th* 
aidpoint of PP* is *xaaln*d: if its sun of equates is lov*r 
than at P», a parabolic interpolation is performed; if 
higher, a point agually distant on th* oth«r aid* of P* 1* 
examined. This process continues until a parabolic 
interpolation can b* performed. Tho us*r is th«n inforaed 
of th* best valu* of R found, together with th* associated 
non-linear paraaetar values, and invited to continue from 
this point, if desired, by again choosinq a search 
technique. If th* us*r is satisfied, however, h* is 
inforaed of the best values found for all paraa*t*rs f and 
the final sua of squares. 

Failure Conditions In NONLXN 

There are two noraal failure conditions wh«ch NONLIN may 
encounter; the usar is inforaed of then by the messages 
SINGULAR MATRIX or STEPS FAILED, and tn* program FLOWNF? 
halts. Singular or near-singular matrices are not uncommon 
in least squares searching. often, however, this failure 
condition has resulted from erroneous network input data 
(see TABLES I, II, and III). The aessaqe STEPS FAILED 
indicates that NONLIN was searching for a minimum in the 
direction of steepest descent (GRAD), but was unsuccessful. 
This failure may result from granularity in the hvpersurface 
which NONLIN is examining. It nay be that the difficulty is 
confined to a small region about the current point P in 
parameter space, and will not recur if the search is resumed 
at a different point. To determine this, after a failure 
message is received, the user must change his trial values 
of the unknown pressures in the flow network. The user is 
on his own if he obtains either of these failure conditions. 

DESCRIPTION OF INTEPPOL 

INTERPOL contains the interpolation and numerical 
differentiation routines used by FLOWNET. Double 3-point 
Lagrangian interpolation is used. In this method, two 
parabolas are fitted to the given data; one to the points 
centered below the argument value and one to the three 
points centered above the argument. A weighted average of 
these two values is the resulting interpolation value. This 
double, 3-point method provides smoother Behavior between 
tabular values than can be obtained with single, 3' point 
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Lagrangian intarpolation. Also, this aathod orovldas a 
function that has • continuous flcat derivative. 

Subprogram 

Following in a liat of subprograms containaa in xhtfppoi, 

A sauna that 

A is an arguaent tabla naaa 
F is a function tabla naaa 
HA is tha langth of A 

X is tha arguaant valua for which tha interpolated 
valua is dssiraa 


Than 

FHTRP (A,F,X,HA) provides tha valua of a raal function, 

DNTHP (A, F, X,HA) providaa tha derivative of a raal function. 

Thasa functions aaka usa of fiva othar functions in 
1MTERP0L. Tha nanas of thasa othar functions ara: fhtfpa, 

FHTFPC, DWTRPC, LIBIT, and TLO. 



APPENDIX B 


SAMPLE PPOBLEM 

An axanpla of tho uaa of tha coaputar prograa PLONMPT 1« 
glvan in this appandix. Xncludad nro tho flow natwork data, 
tho flow path data, tha coapntor torninal output, and tho 
high apood printac output (coaputar output) ralatad to tho 
aanpla problan. Pigura A la a achaaatic raproaantation of 
tha flow natwock uaad for tha aanpla problan. Tha unknown 
praaauraa ara at nodaa 1 through 6 . Tha ayatan inlat 
praaaura la 300 pala and tha gaa uaad la nir. 


25 



Network Data for Sample Problen 


SAMPLE PROBLEM - U.S. CUSTOMARY UNITS 
6PATHSP INLET-8,1, 1,2, 2, 3, 3, 6, 4, 5, 

OUTLET- 1,4, 2 ,6 ,3,9,9,10,5,7, 

FLPTYP® 10* 1, SAME ( 1,2) -4, SAME<2,4)®2, SA MB (3 , 7) *6 , 
NTPATH-10 , SAMEPR-T BEND 

6NODESP NTNODE* 10, FIXBD=6*F, 19*T, HDNET-25*0 . 0, T0=25*70.0, 
PRESS-70.0, 25.0, 15.0, 25.0, 20.0, 15.0, 

14.7, 300.0, 14.7, 14.7 BEND 
6PARAH NPRAT (1 ) -10*8 , DELPMX-0.05, 

PRAT (1, 1) = 1.0, 1.01, 1.05, 1.1, 1.5, 2.0, 3.0, 

4.0 BEND 
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Flow Path Data For Saapla Problaa 


FIRST QUASC SEAL { PATH ONE ) 

BSDATA R1IN-2. 7 , R2IN-0.0, RDIFIN-0.8, WIDTH-0.0, MOLWT-29.0, 
CP-0.24, GAMMA*1. 8 , HU-0.0, SPEED-0.0, CAPV-0.0, XLAB-1.0, 
XTURB-0.25, CONLAM-24.0, CONTRB-0.079, PELAM-2300 . 0, RETUPB-3000. 0, 
PWRSKP-T, PRSSKP-T, NRMSKP-*, PLTSKP»8*T, NOSI-T, SKPH-F, LOSS-0. «, 
IPATH* 1 BEND 
1 

0.01 

FIRST VENT - 4 PIPES , I.D.-0.5 IN., H-0. 25 IN. ( PATH TWO ) 

BSDATA F1IN-1 • 0, R2IN-0.0, RDIFIN-41.0, WIDTH-0.785, MOLWT-29.0, 
CP-0.24 , GAMNA*1.4, HU-0.0, SPEED-0.0, CAPV-0.0, XLAM-1.0, 

XTURB-0. 25, CONLAH-24. 0, CONTRB-0.079, PELAN-2300 . 0, RETUPB-3000. 0, 
PWRSKP-T, PRSSKP-T, NRHSKP-T, PLTSKP-8*T, NOSI-T, SKPH-F, LOSS-0. 8, 
I PATH- 2 BEND 
1 

0.250 

SECOND QUASC SEAL ( PATH THREE ) 

BSDATA R1IN-2.5, R2IN-0.0, PDIFIN-0.8, WIDTH-0.0, MOLWT-29.0, 
CP-0.24, GAMMA-1. 4 , MU-0.0, SPEED-0.0, CAPV-0.0, XLAM-1.0, 

XTORB-O .25, CONLAM-24.0, CONTRB-O. 079, RELAH =2300 . 0 , RETUPB-3000. 0 , 
PHRSKP-T, PRSSKP-T, NRMSKP-T, PLTSKP-8*T, NOSI-T, SKPH-F, LOSS-0.5, 
IPATH-3 BEND 
1 

0.01 

SECOND VENT - 2 PIPES, I.D. *0.5 IN., H-0.25IN. ( PATH FOUR ) 

BSDATA PI IN = 1.0, R2IN-0.0, PDIFIN-41.0, WIDTH-0.785, MOLWT-29.0, 
CP-0.24, GAMMA-1. 4 , MO-O.O, SPEED-0.0, CAPV-0.0, XLAH-1.0, 
XTUFB-0.25, CONLAH-24. 0, CONTRB-0.079, RBLAH-2300 . 0, RETU»B=3000. 0, 
PWRSKP-T, PRSSKP-T, NRHSKP-T, PLTSKP=8*T, NOSI-T, SKPH-F, LOSS-0.5, 
IPATH-4 BEND 
1 

0.250 

THIRD QUASC SEAL ( PATH FIVE ) 

BSDATA R1IN-2.375, R2IN-0.0, PDIFIN-0.8, WIDTH-0.0, HOLWT-29.0, 
CP-0.24, GAMMA-1.4, MU-0.0, SPEED-0.0, CAPV-0.0, XLAM-1.0, 
XTURB-0.25, CONLAM = 24 . 0, CONTRB-0.079, PELAH-2300 . 0 , RETUPB-3000. 0 , 
PWRSKP-T, PRSSKP-T, NRMSKP-T, PLTSKP=8*T, NOSI-T, SKPH-F, LOSS-0.5, 
IPATH-5 BEND 
1 

0.01 
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PIVE EIGHTS PIPE, I. D. *0.625 IN., H-0.313 ( PATH SIX ) 

SSDATA R1IN*1.0, R2IN-0.0, RDIFIN*60.0, If IDTH-O. 982, NOLWT*29.0, 
CP*0.24, GAWHA-1.4, HO*0.0, SPEBD-0.0, CAPV*0.0, XLAH*1.0, 

XTURB-0.25, CONLAH*24. 3 , CONTRB-O. 079, RBLAN-2300 . 0, RFTUPB-3000. 0, 
PNRSKP=T, PRSSKP*T, NRSSKP«T, PLTSKP»8*T, NOSI*T, SKPR*F, LOSS*0.5, 
IPATH*6 SEND 
1 

0.313 

SIX - HALF INCH VENTS AT EXIT, I.D.-0.5 IN., H*0.25 IN. ( PATH SEVEN ) 
SSDATA P1IN=1.0, R2IN-0.0, RDIFIN-24.0, WIDTH*©. 785, ROLWT-29.0, 

CP*0. 24, GAHN A® 1.4, HU*0.0, SPEED*0.0, CAPV*0.0, XLAN=1.0, 

XTURB a 0 . 25 , CONLAM*24.0, CONTRB*0. 079, RELAS*2300 . 0, PFTHPB* 3000. 0, 
PWRSKP-T, PPSSKP*T, NPNSKP*T, PLTSKP*8*T, NOSI*T, SKPH*F, LOSS*0.5, 
IPATH*7 SEND 

INCH AND A HALF LINE, I.D.*1.5 IN., H*0.75 IN. ( PATH EIGHT ) 

SSDATA RlIN=1 . 0 , R2IN*0.0, RDIFIN* 1 85. 67 , WID?H*2.36, HOLWT=29.0, 
CP*0. 24 , GAUM A* 1. 4 , H0*O.O, SPBED*0.0, CAPV*0.0, XLAM*1.0, 

XT0RB*O .25, C0NLAH X 24. 0, CONTPB*0. 079, RELAH*2300 . 0 , RETOPB=3000. 0, 
PWFSKP*T, PRSSKP-T, NRHSKP*T, PLTSKP*8*T, NOSI*T, SKPH=P, LOSS=0.5, 
IPATH=8 SEND 
1 

0.75 

TWO INCH VENT PIPE, I.D. *2 IN., H*1 IN ( PATH NINE ) 

SSDATA R1IN=1. 0, R2IN=0.0, RDIFIN-8.0, WIDTR=3.14, NOLWT®2°.0, 

CP=0. 24 , GAMMA=1. 4 , MU=0.0, SPEED*0.0, CAPV=0.0, XLAM=1.0, 

XTURB-0 .25, CONLAH-24. 0, CONTPB*0. 079, RELAH=2300 .0, RFT0PB=3OOO. 0 , 
PWRSKP=T, PRSSKP=T, NRHSKP=T, PLTSKP=8*T, NOST=T, SKPH=F, LOSS=0.5, 
IPATH*9 SEND 
1 

1.00 

THREE INCH VENT PIPE, I. D. = 3 IN., H=1.5 IN. ( PATH TEN ) 

SSDATA R1IN=1. 0, R2IN=0.0, PDIFIN*304 . 0, WIDTH=4.71, WOLWT=29.0, 

CP=0 . 24 , GAMNA=1. 4 , 110 = 0.0, SPEED=0.0, CAPV=0.0, XLAM=1.0, 

XTORB=0 . 25 , CONLAN=24.0, CONTRB=0. 079, PELAN *2300 . 0 , FFT0PB=3OOO. 0 , 
PWRSKP=T, PRSSKP*T, NPMSKP=T, PLTSKP=8*T, NOSI=T, SKPH=F, LOSS=0.5, 
IPATH = 10 SEND 
1 

1.50 
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Saipl* Problua couputur Ttrninal Output 


OUTPUT PLOW FUNCTIONS FOR THIS ITERATION? 


ENTER (1,2) FOR (CONVERS ATI ORAL, BATCH) 

1 

NEW STARTING VALUES ONLY? 

N 

ENTER NF,NP,NQ,NN, OP QUIT 
SDATA NN 3 6, NP*6 SEND 
ENTER INITIAL NON-LINEAR PARAHRTERS 
SDATA PHI*70.0, 25.0, 15.0, 25.0, 20.0, 15.0 
ENTER (1,2, 3, 4) FOR (GRAD,LIN,SPIRAL,COHP2) 

3 


BmtoDUcroamoFTOE 
ORIGINAL PAGE IS P 


SEND 


R- 8.058E-02 

PHI 3 6.812E 01 2.493E 01 1.485E 01 2.492E 01 1.749E 01 1 

ANSWERS GOOD ENOUGH? 


N 

ENTER (1,2, 3, 4) FOP (GRAD, LIN, SPIRAL, COHP2) 
3 


R“ 3.569E-02 

PHI= 6.172E 01 2.380F 01 1.485E 01 2.045E 01 1.659E 01 1 

ANSWERS GOOD ENOUGH? 

N 

ENTER (1,2, 3, 4) FOR (GRAD, LIN, SPIRAL, COHP2) 

3 


B= 5.086E-03 

PHI 3 6.069E 01 2.279E 01 1.485E 01 2.002E 01 1.594F 01 1 

ANSWERS GOOD ENOUGH? 

N 

ENTER (1,2, 3, 4) FOR (GRAD,LIN, SPIRAL, COMP2) 

3 

F= 1.631E-03 

PHI- 6.093E 01 2.274E 01 1.481E 01 1.936B 01 1.586 01 1. 

ANSWERS GOOD ENOUGH? 


N 

ENTER (1,2, 3, 4) FOR (GRAD, LIN. SPIRAL, COMP2) 

3 

p — 1*4 10E-03 

PHI- 6 1 054E 01 2. 267B 01 1.482E 01 1.914E 01 1.581E 01 1 

ANSWERS GOOD ENOUGH? 

N 

ENTER (1,2, 3, 4) FOR (GRAD,LIN, SPIRAL, COMP2) 

3 


g — 3.o 05 E-04 

PHI- 6^ 004E 01 2.253E 01 1.483E 01 1.907E 01 1.570E 01 1 

ANSWERS GOOD ENOUGH? 

N 

ENTER (1,2, 3, 4) FOR (GPAD,LIN,SPIFAL,COMP2) 


. 5 14E 01 


. 528E 01 


. 529 E 01 


522E 01 


.5188 01 


. 5 1 4E 01 
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3 

R* 1.084B-04 

PHI* 6.013B 01 2.2H8B 01 1.482E 01 1.892E 01 1.8MF 01 1,«HE 01 

ANSWERS GOOD ENOUGH? 

N 

ENTER (1,2, 3,4) FOR (GRAD, LIN, SPIRAL, COHP2) 

3 

R- 6.932E-06 

PHI* 5.991E 01 2.232E 01 1.483B 01 1. 883B 01 1.563E 01 1.«0«E 01 

ANSNERS GOOD ENOUGH? 

N 

ENTER (1,2,3, 4) FOR (GRAD, LIN, SPIRAL, CONP2) 

3 

R* 3.700E-06 

PHI* 5.994E 01 2.232E 01 1.482E 01 1.882E 01 1.563E 01 1. e 0?E 01 

ANSWERS GOOD ENOUGH? 


N 

ENTER (1,2, 3,4) FOR (GRAD, LIN, SPIRAL ,COHP2) 

3 

R* 2.432E-07 

PHI* 5.997E 01 2.233E 01 1.482E 01 1.882E 01 1.563E 01 1.807E 01 

ANSWERS GOOD ENOUGH? 

N 

ENTER (1,2, 3,4) FOR (GR AD, LIN, SPIRAL ,CONP2) 

3 

B* 1.295E-07 

PHI* 5.997E 01 2.233E 01 1.482E 01 1. 882B 01 1.563E 01 1.508E 01 

ANSWERS GOOD ENOUGH? 

N 

ENTER (1,2,3, 4) FOR (GRAD, LIN, SPIRAL, COMP2) 

3 

R* 5.373E-08 

PHI* 5.997E 01 2.233E 01 1.482E 01 1.882E 01 1.563F 01 1.506E 01 

ANSWERS GOOD ENOUGH? 

y 


FINAL RESULTS 

R= 5 . 372823E-08 

PHI* 5. 996709E 01 
2. 232632E 01 
1*482184 E 01 
1.881750E 01 
1.563185E 01 
1 . 506396E 01 

NEW STARTING VALUES ONLY? 

N 

ENTER NF, NP, NQ, NN, OF QUIT 
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BDATJl QUIT-T 6 END 


MASS PLOW BALANCE INFOBNATION AT NODES 1 TO 6 


NODE NO. 

MD(NET) 

MD(MIN) 

MD (NET) /HD (MIN) 

******** 

******* 

******* 


1 

~3. 326E-05 

1.264E-01 

- 2. 632E-04 

2 

-4. 181E-05 

3. 90RE-02 

-1.070E-03 

3 

1 . 992E-05 

8.81 1E-0 3 

2. 260E-01 

4 

5. 662B-06 

5.4441-01 

1 • 040E-05 

5 

-6. 026E-05 

5.444E-0 1 

-1. 107E-04 

6 

-2. 164E-04 

8.732E-02 

-2.478E-03 


ITERATION 1 
NODE NO. DELTAP/P 

1 1. 673E-01 

2 1.198E-01 

3 1 . 202E-02 

4 3.285E-01 

5 2.794E-01 

6 4. 246E-03 
SOLUTION NOT CONVERSED 


BEGIN OUTER ITERATION NO. 2 


OUTPUT PLOW FUNCTIONS POR THIS ITERATION? 

Y 

BEGIN NEXT INNER ITERATION 

ENTER (1,2) FOR (CONVERSATIONAL, BATCH) 

1 

NEW STARTING VALUES ONLY? 

Y 

' ENTER INITIAL NON-LINEAR PARAMETERS 

BDATA PHI=60 * 0, 22.3, 14.8, 18.8, 15.6, 15.0 BEND 
ENTER (1,2, 3, 4) FOR (GRAD, LIN, SPIRAL, COMP2) 

‘ 3 

p -5 1 . 786E“04 

PHI= 5 * 999E 01 2. 230B 01 1.483E 01 1.894E 01 1.563F 01 1.520E 01 

ANSWERS GOOD ENOUGH? 
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N 

ENTER (1, 2 , 3 ,*) FDR ( GR AD, LIN, SPIRAL, CON P2) 

3 

R® 7.277E-05 

PHI- 5.999E 01 2.231E 01 1.481E 01 1.898B 01 1.F62E 01 1.518E 01 

ANSWERS GOOD ENOUGH? 

N 

ENTER (1,2, 3,8) FOR (GRAD, LIN, SPIPAL,COHP2) 

3 

R= 3.787E-05 

PHI- 5.999E 01 2.234E 01 1.481E 01 1.904E 01 1.?63F 01 1.N18B 01 

ANSWERS GOOD ENOUGH? 

N 

ENTER (1,2, 3,4) FDR (GRAD, LIN, SPIRAL, COHP2) 

3 

R~ 8. 3948-06 

PHI® 6 ! 011E 01 2.232E 01 1.481E 01 1.906B 01 1.585E 01 1.815E 01 

ANSWERS GOOD ENOUGH? 

N 

ENTER (1,2, 3,4) FOR (GRAD, LIN, SPIRAL ,COHP2) 

3 

R® 1. 0511-06 

PHI® 6.007E 01 2.234E 01 1.481E 01 1.909E 01 1.566E 01 1.515E 01 

ANSWERS GOOD ENOUGH? 

N 

ENTER (1,2, 3,4) FDR (GRAD, LIN, SPIRAL ,COMP2) 

3 

R= 2.390E-07 

PHI® 6.008E 01 2.234E 01 1.481E 01 1.909E Cl 1.566E 01 1.515B 01 

ANSWERS GOOD ENOUGH? 

N 

ENTER (1,2, 3,4) FOR (GRAD, LIN, SPIRAL, COHP2) 

3 

p~ 8. 751E-08 

PHI® 6l007E 01 2.234E 01 1.4812 01 1.910E 01 1.566E 01 1.515E 01 

ANSWERS GOOD ENOUGH? 

N 

ENTER (1,2, 3,4) FOP (GR AD, LIN, SPIRAL, CONP2) 

3 

R® 1.800E-08 

PHI® 6.008E 01 2.234E 01 1.481E 01 1.910E 01 1.566F 01 1.516E 01 

ANSWERS GOOD ENOUGH? 

t 


FINAL RESULTS 

R® 1 . 800260E-08 


PHI® 6 « 007884E 01 


2.234119E 01 
1 ,4808331! 01 
1.90978 IE 01 
1 . 566009E 01 
1 . 515520E 01 

HEW STARTING VALUES ONLY? 

N 

ENTER NF,NP,NQ,NN*. OR QUIT 
6 DAT A QtJTT=T SEND 


HASS FLO« BALANCE INFORMATION AT NODES 1 TO 6 


NODE NO. 
******** 
1 
2 

3 

4 

5 

6 


BD (NET) 
******* 
-1 . 281E-05 
1. 194E-05 
~5. 729E-06 
1. 318E-04 
1.59 IE-05 
-6.497E-05 


ND (BIN) 
******* 

1.26 IE-01 
3.897E-02 
8.766E-03 
5.445E-01 
5.445E-01 
8.712E-02 


BD (NET) /BD (BIN) 

-1. 016E-04 
3. 064E-04 
-6. 536E-04 
2. 420E-04 
2.923B-05 
-7. 457E-05 


SSSSSSSS»°S» 
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Sample Problem Friiitei Output 


SftHPIH F503IFE - 0. 
PiTH HO. KCDmflH) 


S. C’lCIOKSrT ost? 
SCP3(CU?> PLOT ?F 

*•*•*•*»* ****** 


POCTT. ?»TF" 

• ••• 


1 

2 

3 

f 

b 

7 

a 

5 

ie 


e 

i 

1 

2 
2 

3 

4 
6 


1 

t 

c 


ir 




r 

It 

*> 




n 



J.G22 "iO. FIXED 

******** ***** 


"EEfSUt 1 ' 

******** 


%fT 1»~5 »L«"3 

**»•**»—*** ************* 


1 

F 

P.-0C03 

C? 


02 

2 

F 

''.2300- 

02 

r*~30€€ 

r 2 

i 

F 

J.15CC3 

r i 

r* 

T2 


F 

C.250C' 

C2 


02 

- 

F 

C.23CC? 

or 


02 

c 

F 

0.13CC3 

j* *■* 
.* £. 

*.?C00f 

02 

T 

1 

C.U70E 

C 2 


C2 

6 

T 

C. 37003 

C3 


02 

c 

T 

f-.l4i*3 

0 2 


02 

10 

an 

4. 

0.1*? r -- 

C2 

r m ~ *\ HO to 

02 


C.7C-7" 
t #* * nr 

clcoro 

r^iwsf" 

c.oooo 
r .o-oo r 
r .c*oc 

f 

-.-'r-r 


cc 


O 

C 




65*.?** 

H'riT= 1 C»fc, *e*o 

P?*I= l.C, 1.CC5355, 1.245S95, 7.C55===, •'.s' - , 2.0, 3.0, t.C, 
1.34?5?9, 1. 055559, 1.SC. 2.2, 3.C. #.0, 2*C.O, l.C, 1.CC ce «, 
1.30, 2.C, 3.0, c.c, 2«:.0, 1.;, I.C'-^r’, 1.0a«cc # 1.0=5---9 ( 
L.C, 2*0. C, l.C, 1. 00 c 959, 7.0C95?% 1.0=r5 r 5, 1.:C.. 2.0, 3.*, 
1.C05559, 1.0*5955, 1.C55 e 9 r , 1.5C. 7. % 3.', *.C. 2*C.C, 1.0 
1. C 3 C S59, 1.50, 2.C, 3.C, 2*C.3, 1.0, 1.oc=cta # l.rr^a 

2.0, 3.0, L.0, 2*0. C, l.C, 1.009599, 1.0455'--, 1.C5°^59. 1..-9 
2*C . 0, l.C, 1.005595, 1.0*-=cc, l.ff-tit, 1.30, 2.C. *. r , 


DtiP*X= C.5C2-C1 


2-C.C, 7.0, 7.rror=e 
1.CS trt s, i.ft'cctf 

1.50, 2lc, 3.2 
».£, 2*0. C, l.f 
l.nf'ecj, ( # <i*cc:c 
1 - dctoep ( l.f- 
2.0, 3,0, *.* 
*02*0.'’ 


k?x= e 

60 SL 


BHPRODUCIBILITy OF THE 



QOASI-OIE DIBEISIOIAL COHPBESSIBLE FLO* SMI PSOGFAB 
IAPUT DATA - 


PIES? QOASC SEAL < PATH OWE > 


HI. I ACHES 
2.7000 

F2,I*CHES 

0.0000 

PLOl LEAGTH.IBCHES 
0.8000 

FLOH WIDTH, IACHES 

0.0000 


OOTPOT DATA - 

HI, I ACHES 
2.7000 

■2, IICHES 
3.S000 

no* LP1GTH, LACHES 
0.8000 

FLO* WIDTH, IACBES 
19.8779 

ABBA, II 2 

15.5823 


PO.PSIA 90LIC0LAP WEIGHT ***** 2*** 

20.000 25.000 »***•*•*• 


P3.PSIA 

0.000 


Cr,BT0/lBB-DE3 F 

0.2*0 


A (LABTAAP1 

2*. 000 


PO/P3 

4.000 


SABBA 

1.400 


If.UTIl*) 

i.r» 


TO, DEG F 
70.0 


TISC0SITT,LB-SEC/I*2 

C.0009 


OPPS* LTB*T ?E tl%Bt4I?| 
2300.*' 


LOSS COEF. 
0.50 


SPEED, FPB 
0.00C9 


KITBPBOLEB-7 

0 .'*obp*-01 


T.FT/SEC 

0.9C 


*tT*PBOl' > m 

0.2*00 


10*F* LIBIT *» f?B*FBt*W?| 
3000.0 


PO.PSIA 

300.000 

93,PSIA 

7C.C00 

P0/P3 

4.286 

TO, TEG ? 

70.000 

GAS T3»STA»T,IB-PT/1BB-P 
53.28964 

LOSS COEF. 

0.50 


BOLECOLAE **I3*T 
25.000 

CP, PTB/13B-9F3 5 

0.000 

6ABBA 

1.400 

AISCOSITT, I P-SEC/IE2 
0.2640E-08 


CAtCBlAT* 

••••••••• 


PI"! 



■< 


sprro.FFB 

o.ocoo 

T, FT/SEC 

0.00 


• • 

* * - CBOKIBG FILB THICFAESS • 

* ♦ - TPAWSITIO* PEG TOW • 

* / - T09B0LSAT FLOW * 



QUASI-OBI DIHEBSIOXll COHPFESSIBLE FLOW SEAL P50SFA* 


ftps? ie*'" - a pipes. i.?>.=c.3 tb, h*c.25 tb fr*?a -■*>» 


DATA - 

El, IBCHES 

PO, PSIA 

5°IECUIXF WFIIP* - 

F=«r/S'—* 

1.0000 

2CC.00P 

20. 5f« 


F.2, IBCHES 

P3.PSTA 

CF.FTn/l2S-rt7 F 


O.OOCO 

o.oto 

0.2*0 


FlOS IEBGTH, IBCHES 

FO/«3 

5M»l 

« r-wi»u* 

HI. COCO 

fc.2S€ 

1.CC0 

l.r* 

FlOS HI DTE. IBCHES 

TO, DEG F 

BISfOSTTT, LB-SIT AIK2 

ITBIT *E tISSIB ASf 

0.7650 

70.0 

c.ococ 

23t**.0 


loss COEF. 

SPFEr.*?" 



0.50 

C.0''" r 



f.FI/S'i 

r 


-f-nr.nf r-.-j 

-.252f 


OUTPUT DAT* - 


»r*SFE *wr: f?"=S*?l2tT1 

3 C 7 C.? 


El, IECHES 
1-CCOO 


?C,?SIA 

70-000 


•OIECOIAF F~:"S- 
2P.35E 


>»«« »■ » «« 


B2. IBCHES 
B2.CCC0 


.000 


cp,*Tu/!.?"-rir 7 r 

A .'>«« 


FLOH IEBGTH, IBCHES 
B1.000C 


P0/P3 
2. SCO 


i.aof 


rtoa HIDTH, J5CHES 
0.785C 


TO, DEG F 
70-000 


flSCCSITT.EE-Etr/Tf’ 

r. 2 **o*-''S 


T T*?? 


ASEA.IS2 GAS COFSI A1tr,tB-Pr:/I.EB-F SpEEr. c *"S 

£538. €250 53.2£222 C-CCCC 


E05S COEF. 
0. 50 


*.r:/s r c 

r."r 


•••••••••••••••••••a 

• V 

• • - CF0KI5S PIE* THICKNESS • 

• ♦ - TFAKSITIOB FESICF • 

• / - TBF5012XT FIOH • 


REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 



0«ASI-O»E C I SESSIONAL COFPSESSXELE FlOF 5EAi. PrOGFAH 
IS POT DATA - 


SECOSO O0ASC **•*! t TF3! 


FI, IKCHE5 
2-5000 


PO.»SIA 

70-00C 


•CIFCHI IF ir’ISI* 
2C.f«or 




E2, ITCHES 
C.C030 


P3.PSIA 

o.oco 


CF,BTn/T.3'!-''F7 » 
3.2*0 


latrrtn 

2*.*”*** 


FLOH LEBGTH.rWCHFS 

c.eooo 


PC/P3 

2.8ce 


Sl«* 

l.eno 


1.M 


F1CB HIBTH.ISCHES 
0.0C0C 


TO, DBG F 
“0.0 


*IFCOSTTT,I«-F5C/Tr! 

C.orco 


"5?*- ir*~~ re (u«i«n 

<? 


loss crsr. 
C. 50 


'MPEFO.FTF 

n .C£Ct' 


M*p T »'j! r n 


r.oc 


r |»SFIF!»r» 

0.25"** 


OOIPBT DATA - 


LO»I? U»!T ?F (!B”i;ilT) 

wr.» 


El, I ACHES 
2.5000 


=C,FSIA 

7C.0CC 


*ci*coi*r amir 

2 e .903 


runiiT! 

« >»»« »««« 


E2,X*CHES 

3.3000 

FtOH LF1GTH, IBCEES 
0.8000 

rica *XDTH,IBCHES 
18.2212 

ASEA.X82 

1A.577C 


?3,PSXA 

2 e .COC 

FO/P3 

2.800 

TO,?EG F 
70. COO 

Gf.S COSSTAFT.lB-FT/lr**-' 
53.2890 

10SS COFF. 

C.50 


0. POO 

«I*WA 

1. ***" 

TTsccsrri, is-ssc/t-? 

O.iFAO'-C® 

FrSB%”F 

C.wr« 

l.fF/S** 

f».p« 


* 

• • - CHCrir* r iiF ’•hxcsses* 

• ♦ - r=A*SITXC»* fF^XRW 

* / - '0FE11EFT Fir* 

* 

M»«»M**»»M***»**»M*(*»* 




ffcom we*" - 2 


t .^ o .5 :% 


■= r .25 T* (?*** *** 1*1 


QUASI- Ollc DIHEWSIOBAt C03EF2SS13I E FLOS SEW PFOGFA- 


ISPUT DAI A - 

El.ISCHES 

1 . 00 C 0 


PO.CSIA 

7C.300 


F2,I*CHES 

C.COOC 


P 3 . r SIA 

C.COC 


FiOS 

IEFGTH.HCHES 

*1.0000 

PO/P3 

2.60*1 

FLOS 

MIDI H, I8CHES 

70, DEG T 

•; n ti 


C.“650 

f Ua J 


I.OFF COBF. 
0.50 


•OLECOUF SFT5HI” 

28.?** 

<~ P , rTn / I - , -"' S " - 
C.2*C 

Sl »^» 

i.irr 

fTscc ' in ,' s - 5Er / rt ? 

*V1f«9 

'' I *®* 1 ” 

n .' JOC *' 

»,F?/S*C 

c.or 

10 * A ? 


OUTFOT DATA - 

F 1 .IWCHES 

1.CC0C 

F 2 ,ISCHES 

* 2.0000 

H .08 LE*GTH,I 8 CHES 

* 1.0000 

nos 8 IDTH, IHCHES 

0.7e5C 

AREA.IA 2 

5536. €250 


?C,*STI 

2S . 02C 

P 3 ,*SIA 

15.000 

FO /®3 

1.667 

TO.DEG F 
70.000 

,AF COSSTAIT.IB-F'/IPB-F 
53.35222 

lOS e COFF. 

0.50 


5C1ECBIA5 FFT5HT 

2b. *66 

c '. JW / irv "!? f 

0 . 000 

SAS8A 

1 . tf " 

*IFC5FI*T,L*-F5C/TK2 

**. 26 « 0*‘-»'F 

w ? r53 .? F7 

0 . 0*00 

i .7-/ s*r 

0 . 0 " 


• 

• • - CSOFIF 5 ?IL" T ! IICF 5 EF 5 

• • - TTEAFFITI 3 * «FGI 05 

• / - TDrOOlSST 710* 


*«•••••*» 


»PIWM 

2 *. w ? 

t 

i . r * 

nr* F r iusir*T» 

2300 .-* 

F |*** FF ' , E lr 5*1 

*. 25 CC 

3 * 00.0 


cur ***** 




QtJASI~0*E DIEEFS XOIA1 COSPSESSI3I.E FLOW SFAl PFOGFA" 


TSIF3 QtMSC 'III f *>ATH Fit; I 


IEPOT DATE - 

H1.IKCHES 

2.3750 

E2, INCHES 

C.0000 

Flos iekgth.ieches 

o.eooc 

FLOS aX&TH.IBCHES 
C.CCOC 


OOIFDT DATA - 

F1 # XKCH£S 

2.3750 

F2« UCHES 
3.175C 

FIO 3 LEIG7B, IXCHES 
C.30CC 

flow bidth, isches 

17.*35e 

ABEA, IB2 

13.9*67 


PC.ESX1 POUC81A? FFTII T 

29.C0C 2*.03C 


*«*/**••* 

»»«« «»»»» 


F3.FSIA 

o.ooc 


CP.B'WTTT-EB'S 

C.2*9 


tauitin 

2*.3tSJ 


PC/P3 SAPS A MUril!) 

1.967 1.M*P l. r 5 


TO. DEG * 
79.0 

10S£ COSf. 
C.E5 


FISC0SITT.1P-SEC /’•’I 
C.OOCO 


p.oc"p 


c.oc 


WPT5F LIST" f£ cunn*» 

239r.p 

K|70F*Jrt»* , n 

".7?«r*-«i 


0.2*00 

»0»FF 1X9'? FE 

3O0C." 


PC, PS? A 
25 . 700 

®3,**SIA 

S5.C0C 

P0/F3 

1.667 

TO.PEG F 
70. COO 


GAS C7XS7AS7. 

53.2696* 

lOSS CP'?. 

e.so 


?^l«T8I»5 be; .ft 
29. COO 


r»,?T0/iM-n?5 * 

0.000 

«im 

1.***C 


FI SCOSTTF, 1E-SIC/t?2 
0.26*95-P 6 


SPITD,?** 

C.OOOG 

F.FT/SEC 

O.CO 


• • 


C*1P7I»-ob 


o» 


fact 




* • - CP0KT5G FXia ?HIC**ESS • 

* ♦ - T*»FSI?IO» FFSTO* • 

* / - T0FP1UE5T FlO* • 





QUASI— ONE DIMENSIONAL CONPFESSIELE FLOS SEAL PEOGFAM 


FIVE-EIGHT* CIL D? AI*< 


H=0.?73 fPATH *IT) 


. T .D. 


= 0 . 625 , 


INPUT DATA - 


PI, INCHES 

PO.PSIA 

RC15C3LAF HEIGH"’' 


..0000 

25.000 



R2, INCHES 

?3,FSTA 

CP, F TD/LBN-UES * 

K 1LANINA5) 

0.0000 

0 m f % 0 n 

C.2t" 

2*.7C0 

FLOH LENGTH, INCHES 

FO/P3 

GANNA 

F fLARIBAF} 

60. COCO 

1.667 

1.*C0 

1. r C 

FLOS SIDTH, INCHES 

TO, DEG F 

TISCOSTTT,LT>—ST 0/1*2 

"PPE 3 L’HIA" (?ART*AN) 

0.5520 

7 0.0 

o.ccco 

230C. n 


LOSS COEF. 


K<"U=3U’r»"l 


C .50 




T.FT/S -r N f "M*9DI 3<I7| 

r .oc r -. 25 rir 


LOBE? t’ijtt ;i- 

r 

OUTPUT DATA - 


HI, INCHES 
1.C0PC 


FO.FSIA 
15. COO 


NOLTCTL A = r-IGET 
2S.966 


CALCS’ A""? 

*••••»»»• 


E2, INCHES 
61.0000 


P3, FSI A 

Itt.vnr. 


CP,BTU/I»M-E-*g t 

0.000 


FLOS LENGTH, INCHES 
60. COCO 


PC/P 3 
1.C2P 


SAHRA 

1.C0<? 


FLOS BIDTH, INCHES 
0.5620 


TO,DFG F 
TO.OCO 


f T SC3SrTY,L3-S T, C/'T2 
0 . 26 F 0 E-r c 


AP.EA, IN2 
11666.7186 


GAS CONSTANT, LE-FT/LBR-F 
53.35222 


SPEED, «*• 

r.o r eo 


LOSS COEF. 

0. 50 


T, FT/SFC 
C.0" 






* * - CHOKING FI’R THTCF5FSS • 

» ♦ - Tr ANSI? ION 7FGION • 

* / - TUFEUIEFT ®LON * 


FIOT 


S3 



QUASI-ONE DIMENSIONAL COHPEESSIBLE FLOS SEAL P30GFAH 
INPUT DATA - 


SIX-HALF IECH TrtTS AT *7IT, t.D.=0.5 IK, *=".25 I* fPITB SE»?»| 


FI, INCHES 

PO.FSTA 

HOLECULAP » r I3HT 


1.0000 

15.000 

26.966 

**•*•***« 

F2, INCHES 

P.3,?SIA 

CP,ETU/IM-P53 F 

K a.ART'NAF) 

C.0000 

C.oor 

0.2*0 

2«.0 r 0 

FLOH LENGTH, INCHES 

PO/F3 

GAMA 

s (Lapimhi 

2ft. 0000 

1.020 

i.oru 

*.*c 

FLOH HIDTH, INCHES 

TO, DEG F 

TTSC0STTT,LB-SFC/TN2 

UPr?e UNI* FI f LAMINA 3| 

O.^BSO 

7C.0 

9.0CC0 

7300." 


LOSS f~OEF. 

SPEED, ~®U 

KfM-MLMTI 


' • ' 

U.W'ii) 

0.‘"?' , OF-01 



T,rT/S r C 

iriT^smiLFVl 



O.CC 

r.? c “r 


LONI? ITU'EOLENT) 

3C0C.C 


OUTPUT DATA - 


El, INCHES 
1.00C0 


PG,PSIA HOLECUIAF HPIGH* CAICULATF 

15.000 28. ••**»•••• 


E2, INCHES 
25.00CC 


P3,PSIA CP, BT'VLEH-DEG P 

1ft. 700 0.009 


FLOH LENGTH, INCHES 
24.CCOO 


PO/P3 

1."2C 


GIBS A 
I.eOO 


floh hidth, inches 
0.7e5C 


TO, PEG F 
TO. 000 


fISC0SITT,LB-SEC/IN2 

C.26*9E-C6 


PL O’ 


AREA, IN2 GAS CONSTANT, LB-FT/LBH-E SPEED, ?FN 

1960.3533 53.35222 C.OCOP 


LOSS COEF. 
C .50 


T,*r/s*c 

P.OC 


• * - CHOKING FILE THICKNESS • 

• ♦ - TRANSIT'D* EEGTON * 

» / - TUP FPL EFT FIO* • 

• * 


REPRODUCIBILITY OF THE 
ORIGINAL PAGE IS POOR 



QOASI-ONE DIHENSIONAL COaPPESSIBLE FLOB SEAL PHOGFAH 
INPOT DATA - 


INCH AND A HALF LINE, I.D.=1.5 IN, H=C.7 e IS tPR’'* NtGH T | 


El, INCHES 
1.C0G0 

E2, INCHES 
C.COOO 

FLOB LENGTH, INCHES 
165.670C 


FLOB WIDTH, INCHES 
2. ?60C 


OOTPUI DATA - 

PI, INCHES 
1.COCO 

F.2, INCHES 
106 .£700 

FLOB LENGTH, INCHES 
165. €700 

FLOB BIDTH, INCHES 
2. 3600 

ABEA, IN2 
109467.7500 


PO.PSTA HOtECEtAE BEI3HT »=* ✓*!*•* 

15. COO 28.966 


P3.PSIA 

C.OOC 


CP, ET0/LBN—PE3 * 

0 . 2*0 


Mimw 1 ! 

2*.0C0 


FO/P3 

1.C2C 


GRHHA 

1.COO 


s ataiNANi 

1.CP 


TO, DEG F 
70. C 


nscosiTt,i*-5tc/ir2 

O.OOOC 


P?FE? LIFT- FE ft AFINAS! 
230C.C 


LOSS COEF. 
0.50 


SPEED, FPN 

0.00 f> 0 


KlTBPPOt^NT! 

» ,-»ooCF-'* ' 


f , F-/sr c 
0.00 


C.250C 


LOSE* LIN IT F? f?0=EPLENT| 
■sisp," 


PC.PSIA 

FOLECOLA® NEIGH* 

15. OCC 

2®. 966 

P3,FSIA 

CF.BTD/IBH-DFS 5 

1* .700 

o.coo 

PO/P3 

GAFHA 

1.C20 

1.*00 

TO, DEG F 

NISCOSITT,LB-SEC/IT2 

7C.000 

C. 26*03-0° 


GAS CONSTANT, L3-FT/LBH-E 
53.35222 


SPEED, ®?H 

C.0C00 


LOSS COEF. 
0.50 


», FT/SEC 

0 . 0 '’ 


• 

- * - CHOKING FILN TFICK*FSS 

• ♦ - TEAKS ITIOS EFGIOK 

• / - TBPEOIENT FLOB 

• 

*»»*•*»•*»*»*»»*»»*• **»****»» 


CAICOIATE 

«•*#*•**• 




■*' A 



^ ^ | 

QOASI-ONE DI9ENSIGSAL CCS? 

EESSIELE FLOB SEAL FPOGEAH TBO INCH TEST PIPE, T.D.=2 IN, 

E=1 IN CPATH NI«s) 

INPOT DAT* - 




HI, INCHES 

PO,?SIA 

I'OLECOLAF BETSHT 

F* */**••* 

1.COOC 

15.000 

29.96E 

•*•*••*•* 

B2, INCHES 

?3,HSIA 

CP, BTB/LET-DE” ? 

Kfl*HIN»») 

O.COOO 

C.000 

0.2*0 

2*.0F« 

FLOB LENGTH, INCHES 

PO/P3 

GAMA 

MIWISU) 

8.0000 

1 .020 

i.co 0 

1.00 

FLCS NIDTH, INCHES 

■’’O.DEG F 

NISO'S ITT , L3- SEC/IN2 

0TFF3 L t NTT FE |I»«INAS| 

3.1«C0 

70-0 

0.0000 

23CC." 


LOSS COEF. 


F (T0-301ENT1 


0.50 

O.nOCO 

c.’porr-oi 



T,T r /SEC 

EfMSSOLEN*) 



O.nr 

0.2500 




LONE® II«I? 5E (T0N301EN®) 




370C.0 

OOTPCS DATA - 




HI, INCHES 

PO, *SI A 

HOLECOLAF BFI3ST 

CALC9IATF 

1.0000 

25.0C0 

2P.96F 


E 2, INCHES 

P3.PSIA 

C®,E?0/I9N-’*E3 F 


9.0000 

20. "OC 

0.OC9 

i 





w ELOB LENGTH, INCHES 

P0/P3 

GAMA 


8 .0000 

1.250 

1.*C0 

I 

FLP8 NIDTH, INCHES 

TO, PEG F 

TI SCOSITT , L3-SEC/IF2 

PLOT 1 

3.140C 

70.000 

C.?(CCE-C B 


AREA,IN2 

GAS CONSTANT, L3-ET/13F-E SPEED, 5F8 


251.3274 

53.35222 

0.0000 



LOSS COEF. 

T,*T/S TC 



0.50 

C.O" 






♦ * - 

CHOFING FILH '■HICF5FSS • 



* ♦ - 

TRANSITION FESION • 



* / - 

TTJF BOLEN? FLOB • 








QOASI-OSE DIHEKSIONIL COHPPESSIBLE FLON SEAL PFOGFAF THE EE INCH ¥FN? PIPE, I.D.=3 IN, 
IN POT DATA - 


El, INCHES 
1.COOO 


po,psia sriFcotAS bfish” 

2S.CC0 29.966 


E2, INCHES 
0.0C00 


P3,PSIA 

C.CCC 


CP, BTO/LBN-DEG F 
0.2*3 


FLOB LENGTH, INCHES 
30*. 0000 


PO/P3 

1.250 


GAMA 

1.*00 


FLO* BIDTH, INCHES 
*.7100 


TO, DEG F 
70.0 


TISC0SITT,IB-SEC/IK7 

0.0000 


LOSS COEF. 
0.50 


SPEW’D,? PH 

0.0000 


V,ET/SFC 

C .00 


OOTPOT DATA - 

El, INCHES 
1.0C00 

P2, INCHES 
305.0000 

FLOB LENGTH, INCHES 
30*.C000 

FLCB HIDTH, INCHES 
*.7100 


AEEA, IN2 
2922*3.0375 


PO,PSIA 
20. COO- 
PS, PSIA 
1* .700 

F0/P3 
1.361 

TO. DEG F 
7C.00C 

SAS CONSTANT, LB— FI /LE8-F 
53.35222 

LOSS COEF. 

C. 50 


HCLECDLLR BTISHT 
29.966 


CP, E70/L8N-DEG * 

0.000 

GABHA 

1 .S 0 O 


NISC0SIT7, LE- SEC/IF2 
0.26*09-08 

SPEED, F PH 

C .0000 

N» FT/SEC 
0.0C 


* » - CHOKING FIL» THICKNESS 

* ♦ - TEANSITION EEGION 

* / - T0PE0L5NT ELO» 


H=l. 5 IN f?*TH TEN* 




K (LANINAS) 

2 *. 0 C 0 


1.00 

OPPEF UHTT PS (LAPIFAP) 
23? r .O 


K(T9»B»’! ENT) 

o.-teorr-oi 


^N") 

0.25"0 


LONEF UNIT -s* tTcrpuiEKT) 

V»pr.p 


CALCDIATE 


*L07 


••**•**»**»»»**» *••• **»»*»*** 



ITEFATIOK 2 


PATH *0. 1 P(I*LET)=30C- 

PBATIO HASS FLOW 

1.0CC 0.0000 

1.01C 0.1215 

1.05C 0.2668 

1.100 0.3638 

1.50C 0.6132 

2.000 0.6681 

3.000 0.6736 

4 . 00C C. €707 


PATH HO. 2 P (TWIST! =70.0 


PFATIP 
1.000 
1.010 
1.C5C 
1. 100 
1.500 
2.00C 
3.000 
A. 000 


HASS FLOB 
0.0000 
C. 1099 
C. 2431 
C. 3323 
0.5661 
0.621S 
0.63*6 
0.6344 


PATH HO. 3 P(TWIFT) ="’0.0 


PFATIO 

HASS FLOW 

1.000 

c.oooe 

1.010 

C.2530E-01 

1.05C 

0.5564E-01 

1. 100 

0.7639E-01 

1.500 

0. 1306 

2.000 

0.1439 

3.000 

0. 1473 

4.000 

0. 1474 




PATH TO. * P (I51E‘ r | =2 c .O 


P RATIO 

HASS FLOH 

1.000 

C.CCOO 

1.01C 

C.1680E-01 

1.050 

0.4190E-01 

1.100 

C.5745E-C1 

1.50C 

0. 988SE-01 

2.C0C 

C.10SE 

3.0CC 

0.1128 

4. PCO 

0.1131 


PITH HO. 5 ? (ISLET) =23.0 


St 


PFATIO 
1.00C 
1.010 
1.05C 
1.100 
1.50C 
2.00C 
3. COO 
t.OOC 


BASS FLOK 
C.003C 
0.590CE-C2 
0.1642E-01 
0. 252CE-0 1 
C.4357E-C1 
0.*e3S5-01 
0.500CE-01 
C.5014E-01 


PATH HO. 6 ?CTSLE' r )=15.C 


PFATIO 

HASS FLOH 

1.0CG 

o.ooco 

1.010 

0. 1C5 IE-01 

1.05C 

C. 2343E-C 1 

1.100 

0.3213E-01 

1.500 

0.5533E-C1 

2.00C 

C.S130E-C1 

3.000 

0.6316E-01 

4.CCC 

0. 6330E-C1 




P (OUTLET) =15.0 



15 (OUTLET) =15.0 







PATH SO. 7 


P (ISLET) =15.0 


p(onitrT)=i*.7 


PEATIO 

HASS FLOS 

1.000 

0.0000 

1.01C 

0.3607E-01 

1.050 

0.7952E-G1 

1.100 

0.1065 

1.500 

0.1638 

2.00C 

0.2011 

3.000 

0.2039 

*.000 

C.203* 


PATH SO. S P (ISLET) =15.0 r (OUTLET) =1* . 7 


PEATIO 

HASS FLOS 

1 . 000 

0.9000 

1.010 

C.*e29E-01 

1.050 

C.1Q8* 

1. 100 

0.1*92 

1.500 

C.259S 

2.000 

C. 2902 

3.000 

C. 3015 

t.OGC 

0.3026 


PATH SO. 9 P( ISLET) =25.0 ? (OUTLET) =20. C 


PEATIO 

1.000 

1.010 

1.050 

1.10C 

1.50C 

2.00C 

3.00C 

*.C0C 


BASS FLOS 
0.0000 
C. 1835 
0.395* 

C. 5331 
0.86*6 
0.9056 
0.6378 
0.7629 


PFATIO 

1.000 

1.010 

1.05C 

1.1C0 

1.50C 

2.000 

3.000 

•I.O0C 


(So 


PATH HO. 10 


HASS FLOS 
0.0CC0 
0.279# 
C.61S0 
o.8#ee 

1.##7 

1.5S# 

1.531 

1.632 


P(IKIET)=2P.C 





FIHil FESOLTS 

E= 5.372e23E-C6 

PHI= 5. 9S67C9E 01 
2. 232632E Cl 
1.*6218*E 01 
1. 981750E Cl 
1.563185E 01 
1.506396E Cl 


!WSS F10S 3A1JL5CE ISf CaSiTIO* IT SODES 1 TO £ 

**#»*****»**»*••»*•»*«»•*■*»*»•**•»**•******#** 


BODE SO. SD(KET) EDfTIS) 

«»**•*»* ******* ******* 

1 -3. 326E-C5 -».2feE-C1 

2 -S.1615-C5 3. 5052-02 

3 1. 9S2E-C5 5.E11F-03 

« 5.662E-C6 5.**r4E-01 

* -6. 0262-05 5. ***5-01 

6 -2.16SE-C* 8."’32T-C2 


*1 (TFT) /TOCEir) 


-2.*32E-C« 
-1.<“?CE-C3 
2. 25CE-03 
1.C*C*-05 

-1.107£-r » 
-2.S76E-C’ 


$ 




in r 


ITEFATICH 3 


PITH 10. 1 PfIS!.3?)=300. 


PfOBTlF?) =*«!.«* 


PFATIO 

1.C0C 

1. CSC 
t.-ICf' 
1.500 

2 . COC 
3.000 
• .000 

.003 


HASS FLOH 
C.OOCC 
0. 2666 
C. 3638 
0.6132 
0.6681 
o.e-36 
C. € « 0 f 
0. 67 C7 


PATH SO. 2 PfTSl.rri=cO.O 


?n’m r ’T)=io.p 


PF ATIC 
1.0CC 
1.05C 
1.10C 
1.5CC 
2.000 
3. COC 
3. 16? 

«.occ 


HASS FLOS 

o.coco 

0.2073 

C.263B 

C.4835 

0.5315 

0.5*»35 

0.5*36 

0.5636 


PATH SO. 3 ?fISL 3 T)= 6 C.C r fOtr'lET) = 22. 3 


PES.TIO 

HASS F30H 

1.C0C 

C.OCOO 

1.050 

0.B756F-61 

1.10C 

C.651BE-01 

1.500 

0.1115 

2. COC 

0.1230 

2.666 

C. 1255 

3.000 

C. 1261 

®. COC 

C. 1262 



P1TB SO. * ?(TKL7T»=22.3 


?»RTIC 

HASS FLOW 

1.0P0 

o.oecc 

1.01C 

C. 16*735-01 

1. CSC 

C. 3726E-01 

1.1C0 

C.51105-C1 

1. 082 

C.6738E-01 

1.50C 

0.8607E-C1 

2. COO 

0. 97632-0 1 

3.CC0 

C.1CC7 


PUB SC. e "ITT T.FT)=22.3 


frITIO 

HISS FLOW 

1.000 

0.0000 

1.01C 

C.38S1E— 02 

1.05C 

0. 161C2-0 1 

1. ICC 

C. 22605-01 

1.50C 

C.368CE-C1 

1.5CC 

0. 383GE-C1 

2. CPC 

0.*312F-C1 

3.C0C 

0.»*5?E-01 


PATR *0. € ?<TSI.-?)=U.* 


iunc 

HISS FLOW 

1.00c 

c.ooco 

i.cce 

C. 0*262-02 

1.010 

C. 10385-01 

1.030 

C. 231*2-01 

1.1CG 

0.317*2-01 

1.S0C 

0.5*665-01 

2.00c 

C. 50B6E-C-1 

3.000 

C. 62* 12-01 


si ; > (* ) ?j«f W * p t.l ( ( | f 


?(rortET)=i5. i 


(OP7L2?)=ie.F 


; fOPTI.T?1 = 16.7 


O p 

1c 

5! 

85 



PATH *0 


7 


*<ntaTi = i*.e 


5 (00?1^> = 1*„T 


PPATIO 

SASS F10H 

1.00C 

€.9000 

1.006 

0.32**5-01 

1.C1C 

0. 3563E-01 

1.C5C 

0.7855S-C1 

1.10C 

C. 1072 

1.500 

0.1616 

2.C0C 

0.1586 

3.CC0 

0.201* 


OATH 5C. 6 PITSIST) =15. 1 F fOO-LST) =1*.T 


PFATIG 

BASS F1C« 

1.C0C 

o.oocc 

1.01C 

0. *8505—01 

1.025 

C. 77055-01 

1.05C 

0.1090 

1.1CC 

0.1*55 

1.5CC 

C.2610 

2.00C 

C.291* 

3.0CC 

C. 3C28 



o ^ 

S£ 

eg 


vn 

iw 


PATH 90. 5 ®(ITLiT) =1c.P 


~|G!»TI*-T)=1=.f 


PFATI3 
1.000 
1.C1C 
1.G5C 
1. IOC 
1.2C«. 
1.530 
2.00C 
3.00C 


SASS *108 
C.CCOO 
C.139* 

C. 2575 
C. *012 
C. 5215 
C.6508 
0.6617 
C.6311 



PITH ro. 1C 


? (ISIS':) =15.6 



PFATIC 
1.09C 
1.0 1C 
1.05C 
1.063 
1.1CC 
1.5"C 

2. C0C 

3. C0C 


HASS PIOW 
C.OCOO 
C.215? 
0.67SS 
C.5360 
0.6569 
1.125 
1.262 
1.276 


vn 




mm- W' 




?;kal eesolis 


-= 1.8C02E OE-CS 

fhi= e.oo7eei>E 01 
2. 2341 1?E 01 
1. uSOf 33E Cl 
1 . 90976 IE 01 
1. £660095 01 
1. 51552 OE Cl 


RJ)SS FLOK BALANCE EKFGHKAEEGS *7 NODES 1 TO € 

***************************************** •***» 


POPE SO. 

SD(NET) 

REISES) 

RP»rE7) /BDCEES) 

******** 

******* 

******* 


1 

-1.26IE-05 

1. 261E-C1 

-1.C16E-C4 

2 

1.194E-0E 

3. 697E-02 

3.054F-08 

3 

-5.729S-0f 

S.766E-C? 

-5.536E-CO 

0 

i.3ies-cf 

5.4U5H-n 

2.620' J -C4 

5 

1.5S1E-05 

5.445E-C1 

2. c 23E-f 5 

6 

-6.«i97E-Ce 

5 . "12E-C2 

-7 . 4 e 7*-C5 


£ 




IHE KSTVGKK SOLUTION IS 

###**#*#*>*«** ********** 


PATH SO. 

THHODt 

OOTNCDE 

PIIHIIT) 

P(CnTLIT 1 

llfC OU 

1 

3 

1 

3 on. o 

£C. no 

€.67^- 

2 

1 

V 

ec.05 

19.10 

Q. r fci 6 

3 

1 

2 

£ ".C 8 

22 . ?e 

0 . ‘ 1 2 f ^ 

0 

2 

c 

22.30 

15. 1£ 

0 . E' 12 *-C 1 

5 

2 

3 

22.30 

10 . PI 

0 . "3tC" 5 ' c '-''1 

6 

3 

c 

1*».S1 

10.70 

0.87 6 6 "-02 

n 

3 

c 

1 ».C 1 

10 . "C 

C. ? r 21?- r 1 

8 

c 

1C 

1 = . 1 £ 

10.-0 

C.-"13r-ri 

0 

4 

5 

1 ". 10 

15. ft 

C,;»i 3 

10 

3 

n 

15. fc 

io. ~r 

r. c l'45 


VJl 

Ml 


G £ 

11 

£ es 

r Q 

n3 \5a 


X 
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TABLE I - Variable* in NANELIST/PATKfiV 


Variable 

Hill 

INLET 

OUTLET 

PLOTYP 

NTPATH 

SAKE 

SAMEPB 


nsesriEiisn els 

Intaqer array., Inteqera associated with the inlet 
nodes. *NIFT(,J) is the integer associated with 
th# inlet nod* of th# Jth pa*h. 

Integer array. Integers associated with th» 
outlet n outlet { j) is th# integer associated 
with th# outlet nod# of th# .Tth oath. 

Integer array. Integer# associated with the tvo# 

of flow throaqh each path. If FLOtyp (T)«i, *he 

flow throuqh path T i# d##erib»d bv subroutine 

QUASC. I? FLO"YP (I) *2# then *he flow throuqh ' 

path I is described by a u$#r supplied tlow module. 

Integer. Total number of flow paths in th® 
network. (Idor. + ical r«*h* at a nod® ar« counted 
as one pa^h.) 

Two Hxensional in*<-< 7 e- a-rav. Con*ains *h? (2 

number of equ-' val*n* paths at *he illfi node of 
these rathe. (J,I) «L »“»rs that path number 

I consists of ! equivalent oaths and th a ir inlet 
nods number is J, ( e AFF dees rot hav« to b« 
specified if there are r.o equivalent paths.) 

toqical. If s«t »h°n only *h» nressur® 

ratios for th® firs* flow o**h mus* be input and 
these values a^e used for all *!o» paths in *he 
network. If . FALFE . , *h»n a se* of pressure 

ratios for each flow rath in the network must be inout. 
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:asi2a 

( 39 ) 

n«> 

(9«> 

«,«9) 


REPRODUCIBILITY OF THE 

ORIGINAL PAG* tc 



TABLE II - Variable* in kabbltst/bodesp/ 


vsriabU 

Bus 2£C£tiCti2H 

NTNODS integer. Total rusher of nodes in the network, 

TO Peal array. '"osal temperature “>* the nas a* each 

node. TO (.1) is th* temperature at node .7. 

PXXED Logical array. I * *TXfD(J) ».?Pt7?., the pressur* 
at nod* .1 is specified. If PTXf 5 (.7) * . m*P. , *-he 
pressur* at nods .7 is an unknown. 

HDSBT Peal array. let mass flow it nod# J. T* 10t ?fT(J)*s., 
then ncde .7 is a source or a sink. T* KPeP'tD «o, f 
then no aass is taken from the ne*work a* nod* *7, 

( WDH6 ^ (.7) suet he positive if mass is taken from the 
network at node .7 and neqative if mass is add*d to 
the network at node >7.) 

PRESS Peal array. Contains ♦he values of she pressures 
at each node tr. the network. If there ar° ’?*! 
unknown (non-fired) pressures then the first XT? 
values in the PF^SS array mije£ be the usar supplied, 
initial estimates of these nv unknown pressures. 


Clissslan 


c?m 


(2*) 


(2d) 


(2d) 


ORIGINAL PAGL U 


-r tub 

i * 


* 


\ 
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TtBLF TII - Variables fn wmPLIST/Pmt/ 


Variable 

HIM 

IPPJIT 


PRUT 


DELPHX 


NPX 


rfi££EiDli2C 

Intwgar array. contalrs tha nui*b»r of prassur» 
ratios *o b* usal, a* a givan inlet pressure, to 
qenarata *he mass flow functions for each flow oa*h. 
If NPP>T(!)»w, than tba de*ere*e mass flow function 
for path I is generated from w values of the 
pressure ratio. 

Two dimensional r«al array. contains the values 
of the pressure ratios tha* will he used *o 
ealsulnts the mass flow functions for each flow 
path. p»XT(L,. 1) is th» Ith vain* of tha pressur* 
ratio for tha Jth path. For each flow path, tha 
first valu* In tha pressur* ratio array must ha 
1.0, and tha oressure ratio valuas oust ba 
monotonia increasinq. 

Used as a convergence criterion. IT th* seal*'! 
pressua at aach non-fixad node changa* by lass 
than tha valua of DELP5T from ona outar 
itaration to tha next than tha solution is eon- 
sidared converaed. 

Integer. The number of unknown orassuras in tha 
flow network. (This is calculatad by F10VK*?) 
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Ulsinsisa 

(?«) 


( 10 , 50 ) 


REPRODUCIBILITY OP THE 
ORIGINAL PAGE IS POOR 



PMMKP Logical variable. V 

calculations Involving poser 
art skipped 

NIIKSKP logical variable. T? .Ten* 1 ,, 

ncrmalintd values o* sealing dam 
fore# and canter of pressure are 
Skipped. 

PRSSKP logical variable, T? sat .TP'tF. 

printout o* distribution*, acros* 
fact of «»al vi 11 b» o »it*»d. 

PLTSKP Array of tight logical variables. 

If all etgh* art s*t .’’•"its. no 
plot* art mala. 

SKPH Logical variablt. If ,‘*PlfE. no 

film thickness data will b* teal. 

NOSI logical variable. If pet . T°tJF. , 

. input, internal calculation* and 

output ar* in »».S. units. I f e*t 
•FALSE. , input, internal 
calculation*, and output are in sr 
units. 

LOSS Entrance velocity loss coefficient 

IPATH Intagar assiqned to the fiov path. 

(Usad only to identify *he data) 


REPRO! I’CHJLITY OF THE! 
ORIGIN A., *-AGE IS POOR 


.T*g*. 

.tpilf. 
a* .TPitp. 
. FALSE. 
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wits If - Fariablem In w»H*LXWT/8B«W/ 


variable 

HUS 

Mil 

P2IH 

RDIPIM 

WIDTH 

HOIWT 

CP 

HO 

a Min a 
SPSS!' 
CAPV 
XLAn 

COKLAM 

XTUFB 


CONTHB 


RELAM 

RETURB 


Iilas S2 H22 

SSSSZiStiSH I2Z 1121 12U2Z1 

Inntr raSlua of seal 
Outar radio* of s*al 
Flow length 
H«an flow wM*h 
Nolacular weight 
Spacific heat 

Reservoir viacositv. "ha program 
will calculate M'l *or air hu* not 
for other game? 

Patio of so*eiflc h*»*« 

eotaticnnl velocity 

Peal fac» sceed 

Exponent in friction facto’- - 
Reynold* number relation for 
laminar flow (eg. in of ref 1) 

Conatar 4, in *ric*ior. ^actor - 
Reynold? number relation for 
laminar flow 

Exponent In friction factor - 
Reynolds number relation for 
turbulent flow 

Constant in frfe-ion. r acto T " ■ 

Reynolds number relation 'or 

turbulent flow 

Maximum esvnolrts number for 
laminar flow 

Minimum "eynolds number for 
turbulent flow 
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ORIGINAL PAGE IS PQOf 



t»pl F V - P*ti 'or 


*»»♦*©*** o' Hour* * 


Mnl*« 

*nET< 3 ) 
OIIT!,^ H) ■*# 

tio*yo{t) si 


mTH8P 0 BTT-M 2 J-J. 

*Lf>»v» (i )" 1 , n HI* 1 # 

5 *»E< 1 # 2 > ■? ' 

6N0DCSP H^HODE-tt, 

mw.ms*., v.T^.. 

HppjtTO*?, SIME 15 **.* 0 ^ * » . . n 

mni.i, m, i.S >.». 


1U. ^ 


ttkfiM 


REPRODUCIBILITY OP THE 
ORIGINAL PAGE IS POOR 


r.*n»> 

'.EWD 
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Sample Flow Network 
FIGURE 1 




Read Initial Estimates 
of the Unknown Pressures and 
All Network Information 


■ 4 I OUTE R" I OUTE R» 1 



ale r 

Call QUASC to 
calculate mass 

Yes 


No^ 

flow function 



for this Path 





L 


. Path No. 

W Path No. +1 


I 


Stop 


[Ko 

Is this the Last 
Flow Path in the 
Network? 

Yes 

if I INNER- 1 


Use Least-Squares 
Search Technique to 
Minimize the Sum of 
Squares Residuals 

if 


Is the Sum of 
Squares Residuals 
Minimized? 

Yes 

5' 

No 

Are Values of New 
Pressures Nearly 


Rquui vu Values oi 

Old Pressures? 

^ Yes 


Output Network 
Solution 


Call ADDFIO to) 
Calculate Massi 
Flow Function ; 
for this Path ! 
1 4 


I INNER" I INNER* 1 

w 

No 


REPRODUCIBlLITy OF THF 

ORIGIN AT- PAO" "* 


Functional Flow Chart for the Program Flownet 


FIGURE 2 














Flow Network for Saaple Proble» 


FIGURE 4 






